Solutions to the restricted three-body problem with variable mass 
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We look for particular solutions to the restricted three-body problem where the bodies are allowed 
to either lose or gain mass to or from a static atmosphere. In the case that all the masses are 
proportional to the same function of time, we find analogous solution to the five stationary solutions 
7-H ■ of the usual restricted problem of constant masses: the three collinear and the two triangular 

t-H ' solutions, but now the relative distance of the bodies changes with time at the same rate. Under 

, some restrictions, there are also coplanar, infinitely remote and ring solutions. 
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I. INTRODUCTION 



The study of the two-body problem of variable masses began practically with the work of Gylden [8|, who wrote 

! , the differential equations of motion for the problem. The first integrable case to Gylden's equation was give n by 

Mestschersky fl5j . for a specific mass variation law. This mass variation law, and its following generalization |16j . are 
Q-ij known as Mestschersky laws. After Mestschersky's contribution, the physical meaning of the problem became clear 
^ fin and it is known as Gylden- Mestschersky problem. 

Jeans Q, by studying the orbits of binary stars, found a m ore general mass variation law that was based on the 
relation between mass and luminosity of the stars presented bv lEddingtonl in the same year. Mestschersky's laws are 
G special cases of Jeans' law. 

Gelfgat [|| considered a different mass variation law. Q investigated the problem using a differential equation 
I ■ transformation method. Also a number of approximate analytic solution were found, e.g., [17| . [l2j studied the 
, particular problem where the total mass is constant, which can be applied to conservative mass transfer in close 
fT) ■ binary systems fl3j |. 

0^ | The Gylden-Mestschersky problem can also be generalized to include the restricted three-body problem. In this 
approach, it is assumed that the two heavier bodies have their motion determined by the Gylden-Mestschersky 
equations. Thus, one have to deal only with the motion of the third body, which does not affect the main bodies 
motion. It was shown Q that this problem presents particular solutions that are analogous to the stationary solutions 
of the classical problem of constant masses: the three collinear solutions L-\ to L 3 and the two triangular solutions 
L4 and L5. Many other particular solutions have also been found [U I2I Hol. HT1| . S ince then, further characteristics of 

T 7"! ■ this problem have been studied, for example, Luk'yanov Singh and Leke |l9j . 

J> \ Besides the Gylden-Mestschersky problem, there are many different cases of two-body problems with variable mass 
[T§ |. These can be classified according to the presence or not of reactive forces, to the variation of the mass of just one 
or both of the bodies, to whether the bodies move in an inertial frame or not and so on. In this paper, we consider 
the specific case where the three bodies move in an inertial frame, within a static atmosphere, from which they absorb 
mass or to which they lose mass. We search for the five particular solutions analogous to the classic case ones, L\ to 
/..-,- 

In the restricted three-body problem, the motion of the two primary body is determined a priori. In the papers 
mentioned above, this motion is determined by the Gylden-Mestschersky problem [l6j], whereas in the present paper, 
the primary bodies move in a static atmosphere, from which they absorb or lose mass. Therefore, reactive forces are 
present. Fact not considered in the Gylden-Mestschersky problem. 
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II. PARTICULAR SOLUTIONS 



The equation of motion for a body whose mass depend on time, m (t) , see for instance Sommerfeld [20( , is 

F = mv + (v — u) m, (1) 

where F is the sum of all the forces acting on it and v is its velocity, both measured in an inertial coordinate system. 
Also, u is the velocity of the center of mass of the absorbed mass immediately before its union with the body (or of 
the ejected mass immediately after its ejection). The overdot denotes, as usual, derivation with respect to the time 
variable. 

In the present context, there are two special cases of equation (TTJ) to be considered. The first one is when the mass 
is ejected with the same velocity of the body at any moment (u = v), i.e. mass ejection does not produce reactive 
forces. This case can be used to study the motion of a body ejecting mass isotropically (or radiating energy, since the 
total reactive momentum would be zero). In this case Eq. (JlJ takes the traditional form F = mv. This case includes 
the Gylden-Mestschersky problem. 

The second case, which we considered in the present paper, is the one where u = 0, i.e. the particles are at rest in 
an inertial coordinate system. This case can be used to study the orbits of a star moving through a static atmosphere, 
whose particles attach or detach to the star as it moves. In this case, equation ([1]) reduces to 

F = 77iv + 777V = — (mv) . (2) 
at 

According to ©, the equation of motion for a problem of two bodies of varying masses mi (t) and mi (t), exchanging 
mass with a static atmosphere surrounding them, is 

d . . . Gm\m 2 

dt (miri) = 7^^ 

d Gmi_m 2 

dt (m2r2) = ^^ r ' (3) 

where ri and r 2 are the position vectors of the two bodies, r = ri — r 2 , r = |r| and G is the gravitational constant. 
We note that the atmosphere around the bodies does not cause any drag forces on them. The only appreciable effect 
of this atmosphere is to work as a source or sink of mass for the bodies. From Eqs. ([3]) we get, 

^- (mir i + m 2 r 2 ) = 0. (4) 
dt 

The quantity in parenthesis can be set equal to zero in an appropriated inertial frame. It follows, ri = (rri2/M)r, 
where M = mi + 7772, can be used to cast equation (J3|) as 

where fi(t) — m\m 2 jM is the reduced mass of the problem. In the classification table of the different two-body 
problems with variable masses by Razbitnaya [l8j |. the present physical problem is listed as number 3. This problem 
also coincides mathematically with problems 14 and 15 (the later, known as Gelfgat-Omarov problem). Eq. ([5]) can 
be decomposed, in polar coordinates (r, 9), in the relations: 

d k 2 GMfi 

-77 (W = — — > (6) 

at [j,r a r z 

fir 2 e = k. (7) 

The coordinate r is the distance between the primary bodies, 9 is the angle between the line passing through their 
centers and a fixed line in the plane of motion, and 9 is the corresponding angular velocity. The constant k is the total 
angular momentum of the system. Equation together with the definition of the center of mass position vector, 
777,2 r 2 , also yields the relation for the time dependence of R cm 

• d /mi \ , , 

R - -jt(w) v (8) 



which does not perform, in general, an inertial motion. 



3 



We shall restrict ourselves to the case where the masses m\ and m 2 vary arbitrarily, but their ratio remains constant. 
In this case the center of mass of the system moves inertially. Following Luk'yanov [ll[ , the time dependence of the 
masses is described by a positive function u (t) : 

mi = toiom (t) , m 2 = m 2 gu (t) , M — Mqu (t) . (9) 

Furthermore, we have: 

(J, = (J, U (t) . (10) 

The constants m 10 , m 2 Q, Mq and are all positive, with M Q — m 10 + m 2 o and /j,q — m w m 2 ^/MQ. 

Now we introduce a barycentric rectangular coordinate system that rotates comoving with the main bodies, with 
angular velocity u) — uz. The main bodies move in the xy plane, in such a way the x axis always passes through the 
center of the main bodies, at x = x\ and x = x 2 . 

In such a coordinate system the equations of motion for the third (lighter) body are, 

— (mx) — F x + muy + 2muy + muy + mui 2 x, 

d 2 

— (my) = F y — mujx — 2muix — miux + mui y, 

j t {mz) = F z , (11) 

where m (t) is the third body mass and F x , F y and F z are the x, y and z components of the gravitational forces F, 
due to m\ and m 2 . It is instructive to c ompare these two last e quat ions with the corr espo nding equations for the 
motion of the third body in references like iBekovl ()l988l and ll99ll ) and lLuk'vanovl (|l989al and ll989o ). Equations (fTTj) 
have the additional terms muiy and — mux. These terms arise from the fact that the motion of the third body obeys 
equation ([5]). Equations (ITT]) also appear in Bekov et al. 0. 

Following Luk'yanov [11| . we first write the distance between the primary bodies as r — roR(t), where R (t) is a 
positive function, so we have to = ujq/uR 2 (equations ([7]) and (|LT|) '). Then, we define the mass parameter v — TO2o/-^o 
(0 < v < 1/2, without losing generality), and choose units such that ro = 1, Mq = 1 and ljq = 1. Therefore, 

x\ = —vR (t) , 
x 2 = {l-v)R(t). 

Furthermore, we shall assume that the mass of the third body also varies as m = m$u (t), where mo is a positive 
constant. In the new units described above, the equations of motion (fTTj) reduce to, 



dt 



(ux) = — (1 — v) (x + vR) 

r i 

Gu 2 

-)±— v (a, _ (i- v ) R ) + 
2y 2yR x 



R 2 R 3 uR 4 ' 

d . Gu 2 „ . Gu 2 
dt {uy) = -— {l - v)y -^f 

2x 2yR y 

d , , \ Gu 2 , Gu 2 

— (uz) = ; 

dt r 



— {uz) = —(l-v)z —vz, (12) 



and 



where n = \J (x + vR) 2 + y 2 + z 2 

r 2 = \J {x — (1 — v) R) 2 + y 2 + z 2 . The constant mo cancels and it does not appear in equations (|T2"j) . 

We look for stationary solutions to equations (fTS]) of the form x — £R(t), y = r]R(t) and z = (R{t), where £, 77 
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and £ are constants. By substituting these three relations for x, y and z into equations (|12[) . we get 



'df 



'df 



dt 



uR) = - 



uR 



p\ 



V 



iR 3 ' 



Gu 2 (1 + i/) 

i? 2 p? 
Gm 2 ^(C + ^- 1) M J_ 

" ^ 2 Pi 

Gm 2 (1 - i/) 77 
"R 2 

Gu 2 vn 



ui? = - 



Cm 2 (1 - v) C _ Gu 2 K 
# 2 i? 2 p 3 ' 



(13) 



where p x = y(£ + ^) 2 + r/ 2 + C 2 an d 



P2 = y (£ + v — 1)' + 77 2 + C 2 - The relative motion of mi and W2 is given by i? (t). For this class of solutions the 

geometrical configuration of the system is similar at any moment. In order to find the time derivative of uR, Eq. © 
is written in the new units as, 



1 



uR 3 R 2 ' 

Finally, we obtain the form of equations (|13[) that allows us to find £, rj and £ 



(14) 



(1 -*)(£ + !/) 1/(^ + ^-1) 



uR 3 



Pi 



Pi 



77 1- 



1 - 1/ 



1 - GRu A 1 



P? 



Pi 



= o, 
= o, 

= 0. 



(15) 



To find the first two equations, we have divided by the nonzero quantity ^jr- If we take £ = (solutions in the 
plane of motion of the primaries), the first two equations are exactly the corresponding equations for the restricted 
three-body problem of constant masses. This means that all the stationary solutions of the constant masses problem 
are also present in the problem of variable masses discussed here, namely the three collinear solutions L\ to L3 (the 
three masses aligned) and the two triangular solutions L4 and L5 (the three masses at the vertices of an equilateral 
triangle), but now the relative distance of the bodies change with time at the same rate. 
In order to look for coplanar solutions, we take rj = 0. Then, equations (|15|) reduce to: 



5 .3 -3 U ' 



P\ 



GRu 1 



Pi 



Pi 



n 3 



(16) 



The existence of coplanar stationary solutions in the equations above requires 



GRu 3 



constant = n. 



(17) 



In other words, the primar y b odies must perform a particular motion, determined by a particular time variation rate 
for the masses masses. In a similar restriction is found. Whereas in that work the restriction leads the masses 
to vary according to Mestschersky unified law, here, in general, one must solve numerically equations (j!7[) and (| 14[) . 
From equation (fTo) we have that the possible values for k are restricted by the relation k > 1. 

Once equation (fi~7|) holds (with k > 1), equations (fl"6)) determine the existence of coplanar solutions [Toj]: L 6 and 
L-j, and also L 8 to Ln, when the value of k is less or equal to a certain value that depends on v. In the limiting case 
k — > 1, coplanar solutions do not exist, but the infinitely remote solutions L± oc appear, with £ = 7/ = and £ = ±00. 
Solutions L\ to L5 are still present. 
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III. RING SOLUTIONS 



From our previous relations it is straightforward to verify the existence of ring (Lq) solutions. We consider now the 
collinear three-body problem with variable masses 0. This is the particular case, where the primary bodies move in 
a straight line passing through their centers. The motion of the primaries is still given by equation ([5]) and the ratio 
of their masses is still constant in time. We can obtain the radial equation of motion by setting k = in equations 
© and 0: 

with 9 = constant. 

Since there is no need for the rotating frame of reference used before, we obtain the equations of motion for the 
third body by setting u) = in equations , thus returning to the inertial frame. Hence, the analogous of equations 
(fUZj) and (flU|) will present only the first two terms of the right-hand side of each of them. Furthermore, we notice, by 
symmetry, that it is enough to restrict ourselves to the plane z — and rotate the found solutions around the s-axis 
for all the spatial solutions. In the units used before, equation (fl"8|) is cast as: 

d ( ■ \ Gu 2 , „. 

jt H = -I? (19) 



By substituting the above equation in the analogous of equation (fl~3|) . we get exactly the first two of equations (|15l) . 
Therefore, we confirm the existence of solutions Li to L% and also solutions L Q (a ring around the x-axis, generated 
by rotation of solutions and L5). 



IV. CONCLUSIONS 



The presence of a static atmosphere (i.e., a source or sink of mass) in this three-body problem allows the center of 
mass and the relative motions of two main bodies to be different from that ones of the Gylden-Mestschersky problem. 
These motions can only be determined if the time variation of the masses is known. 

In summary, we found that our variable mass three body problem has also the five Lagrange solutions of the classic 
three-body problem, L\ to L5, but now the relative distance of the bodies change with time at the same rate. In the 
particular case in which the problem is collinear, triangular solutions L4 and L5 produce a ring solution Lq- 

The sufficient condition for the existence of these solutions is that the ratio of their masses be constant in time. 
Obviously, this restriction excludes many physically reasonable models for the mass variation, e.g., the ones where 
the rate of change of the masses depends on the masses themselves. 

The coplanar solutions L e to in and the infinitely remote solutions L± oa , that appear when the motion of the 
primaries is given by the Gylden-Mestschersky problem, are also present, under the very same additional conditions 
on the parameter k found in that case. 
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